Temperature-measuring method and distributed optical fiber temperature sensor

ABSTRACT

A temperature-measuring method comprises inputing laser pulses into an optical fiber to be measured and measuring temperature distribution in the fiber from the ratio of the amplitudes and the delay time of Stokes light and anti-Stokes light contained in the return beam from the optical fiber, wherein the temperature distribution is measured by using equations: ##EQU1## where T(x) is temperature to be measured, Θ is the reference temperature, R&#39; (T) is the relative ratio of amplitude at the measuring point, R&#39; (Θ) is the relative ratio of amplitude at the reference temperature point, k is Boltzmann&#39;s constant, h is Planck&#39;s constant, c is velocity of light, ν is Raman shift, α is distance, wherein the attenuation difference α is represented by a function α{T(τ)} dependent on temperature T(τ) at the measuring point τ, and the exponential function in which the exponential part is represented by a value obtained by integrating the function α{T(τ)} with the variation dτ in distance, is employed as a correction factor.

This is a division of application Ser. No. 534,233, filed on June 7, 1990, now U.S. Pat. No. 5,054,935.

BACKGROUND OF THE INVENTION

The present invention relates to a temperature measuring method employing a distributed optical fiber temperature sensor. More particularly, it relates to a distributed optical fiber temperature sensor and a temperature-measuring method, whereby the influence of the attenuation difference in the optical fiber to be measured, which creates a problem when the ratio of the Stokes light and the anti-Stokes light is taken, is corrected so that the temperature distribution in a long distance optical fiber can be measured with low precision.

A block diagram of a conventional distributed optical fiber temperature sensor is shown in FIG. 5. Laser pulses generated from a laser pulser 10 of a light source, are input into an optical fiber 12 to be measured, and Raman back scattered light generated in the optical fiber 12 is returned to the incident end. Such Raman back scattered light is introduced by an optical directional coupler 11 to a measuring apparatus, whereby the Stokes light and the anti-Stokes light in the Raman back scattered light are separated by a filter 13, detected, and converted, respectively, to electrical signals in proportion to their associated amplitudes by the respective photo-electric converters 14 and 14'. Such electrical signals are amplified by the respective preamplifiers 15 and 15' and subjected to a prescribed number of averaging treatment operations by an averager 16. The resulting average-treated signal is transmitted to a signal processing unit 17, whereby the delay time and the ratio of the amplitudes of the Stokes light and the anti-Stokes light are calculated, and the temperature distribution is output.

In the calculation of the ratio of the amplitudes of the Stokes light and the anti-Stokes light, it is difficult to obtain the ratio of the absolute amplitudes. Therefore the ratio of the amplitudes of the strokes and antistrokes light is obtained from the ratio of the relative amplitudes (hereinafter sometimes referred to as "the relative ratio of amplitudes") by the following equation by Dakin (Dakin, J. P. Pratt, D. J., Bibby, G. W., Ross, J. N. "Temperature distribution measurement using Raman ratio thermometry", SPIE Vol 566, Fiber Optic and Laser Sensors III 249 (1985)): ##EQU2## In the above equation, T is the temperature to be measured, Θ is the reference temperature, at a reference measuring point R'(T) is the relative ratio of amplitude at the measuring point, R'(Θ) is the relative ratio of amplitudes at the reference temperature point, k is the Boltzmann's constant, h is the Planck's constant, c is the velocity of light, and ν is the Raman shift.

The incident laser pulses from the light source, such as a semiconductor laser, naturally undergo attenuation during transmission through the optical fiber to be measured. The influence of this attenuation can be avoided by taking the ratio of the Stokes light and the anti-Stokes light. However, the influence of the attenuation difference between the Stokes light and the anti-Stokes light, when the Raman scattering is generated by the incident wave and the scattered light is transmitted backwards, can not be avoided simply by taking the ratio. Therefore, a correction is required which takes the influence of this attenuation difference into consideration. If such a correction is not conducted, even if a uniform temperature distribution is measured, a substantially inclined temperature distribution will be output as shown in FIG. 4. Further, the attenuation difference itself has a temperature dependency, and the influence of the temperature dependency will be substantial in the case of a measurement taken over a long distance. As a result, accurate measurement with high precision has been difficult to obtain.

SUMMARY OF THE INVENTION

The present invention has been made to solve the above mentioned problems. According to a first embodiment, the present invention provides a temperature-measuring method which comprises inputting laser pulses into an optical fiber to be measured and measuring the temperature distribution in the fiber from the ratio of the amplitudes and the delay time of the Stokes light and the anti-Stokes light contained in a return beam emitted from the optical fiber. The temperature distribution is measured by using an equation: ##EQU3## wherein T(x) is the temperature to be measured, Θ is the reference temperature, R'(T) is the relative ratio of amplitudes at the measuring point, R'(Θ) is the relative ratio of amplitudes at the reference temperature point, k is the Boltzmann's constant, h is the Planck's constant, c is the velocity of light, ν is the Raman shift, α is the attenuation difference in the optical fiber between the Stokes light and the anti-Stokes light, and x is the distance. In the equation the attenuation difference α is represented by a function α{T(τ)} which is dependent on the temperature T(τ) at the measuring point τ, and the exponential function, in which the exponential part is represented by a value obtained by integrating the function α{T(τ)} with the variation dτ in distance, is employed as a correction factor.

According to a second embodiment, the present invention provides a temperature-measuring method which comprises inputting laser pulses into an optical fiber to be measured, and measuring the temperature distribution in the optical fiber from the ratio of the amplitudes and the delay time of the Stokes light and the anti-Stokes light contained is the return beam from the optical fiber. The temperature distribution is measured by using first and second correction equations as set forth below

First correction equation ##EQU4## In the first correction equation T(x) is the temperature to be measured, Θ is the reference temperature, R'(T) is the relative ratio of amplitudes at the measuring point, R'(Θ) is the relative ratio of amplitudes at the reference temperature point, k is the Boltzmann's constant, h is the Planck's constant, c is the velocity of light, ν is the Raman shift, α is the attenuation difference in the optical fiber between the Stokes light and the anti-Stokes light, x is the distance, and u is a real number of at least 0, wherein the attenuation difference α is represented by a constant value, and the exponential function, in which the exponential part is the product of the attenuation difference α and the distance x, is employed as a correction factor.

Second correction equation. ##EQU5## In the second correction equation T(x), Θ, R'(T), R'(Θ), k, h, c, ν, α and x are as defined above, with respect to the first correction equation wherein the attenuation difference α is represented by a function α{T(τ)} which is dependent on the temperature T(τ) at the measuring point τ, and the exponential function, in which the exponential part is represented by a value obtained by integrating the function α{T(τ)} with the variation dτ in distance, is employed as a correction factor; In the second embodiment, the temperature distribution [T_(n) ^(o) ], (wherein n is an integer of at least 1) obtained by the first correction equation, is substituted into the exponential part of the second correction equation, to obtain temperature distribution [T_(n) ¹ ]. Subsequently, temperature distribution [T_(n) ^(n) ], (wherein r is an integer of at least 1) obtained immediately before, is successively substituted into the second correction equation, to obtain temperature distribution [T_(n) ^(r+1) ] at r+1 times. As a result, the temperature distribution [T_(n) ^(r) ], obtained when the error ##EQU6## measuring points. becomes not higher than a predetermined value at r times, is taken as the measured value.

According to a third embodiment, the present invention provides a distributed optical fiber temperature sensor which comprises a light source which inputs laser pulses into an optical fiber to be measured, an optical directional coupler which introduces the return beam from the optical fiber to a signal processing unit, and a signal processing unit which detects the Stokes light and the anti-Stokes light contained in the return beam and which measures the temperature distribution in the optical fiber from the ratio of the amplitudes and the delay time of the Stokes light and the anti-Stokes light. In this third embodiment, the signal processing unit measures the temperature distribution by using an equation: ##EQU7## wherein T(x) is the temperature to be measured, Θ is the reference temperature, R'(T) is the relative ratio of amplitudes at the measuring point, R'(Θ) is the relative ratio of amplitudes at the reference temperature point, k is the Boltzmann's constant, h is the Planck's constant, c is the velocity of light, ν is the Raman shift, α is the attenuation difference in the optical fiber between the Stokes light and the anti-Stokes light, and x is the distance.

According to a fourth embodiment, the present invention provides a distributed optical fiber temperature sensor which comprises a light source which inputs laser pulses into an optical fiber to be measured, an optical directional coupler which introduces the return beam from the optical fiber to a signal processing unit, and a signal processing unit which detects the Stokes light and the anti-Stokes light contained in the return beam and which measures the temperature distribution in the optical fiber from the ratio of the amplitudes and the delay time of the Stokes light and the anti-Stokes light, for this fourth embodiment. The signal processing unit measures the temperature distribution by using first and second correction equations as set forth below.

First correction equation: ##EQU8## In this first correction equation T(x) is the temperature to be measured, Θ is the reference temperature, R'(T) is the relative ratio of amplitudes at the measuring point, R'(Θ) is the relative ratio of amplitudes at the reference temperature point, k is the Boltzmann's constant, h is the Planck's constant, c is the velocity of light, ν is the Raman shift, α is the attenuation difference in the optical fiber between the Stokes light and the anti-Stokes light, x is the distance, and u is a real number of at least 0, wherein the attenuation difference α is represented by a constant value, and the exponential function, in which the exponential part is the product of the attenuation difference α and the distance x, is employed as a correction factor; and

second correction equation: ##EQU9## In the second correction equation R(x), Θ, R'(T), R'(Θ), k, h, c, ν, α and x are as defined above, with respect to the first correction equation the attenuation difference α is represented by a function α{T(τ)} which is dependent on the temperature T(τ) at the measuring point τ, and the exponential function, in which the exponential part is represented by a value obtained by integrating the function α{T(τ)} with the variation dτ in distance, is employed as a correction factor, in this fourth embodiment. The temperature distribution [T_(n) ^(o) ] (wherein n is an integer of at least 1), obtained by the first correction equation, is substituted into the exponential part of the second correction equation, to obtain temperature distribution [T_(n) ¹ ], subsequently, temperature distribution [T_(n) ^(r) ] (wherein r is an integer of at least 1) obtained immediately before is successively substituted into the second correction equation, to obtain temperature distribution [T_(n) ^(r+1) ] at r+1 times, and temperature distribution [T_(n) ^(r) ] obtained, when the error ##EQU10## (wherein m iS the number of the measuring points), becomes not higher than a predetermined value at r times, is taken as the measured value.

BRIEF DESCRIPTION OF THE DRAWINGS

In the accompanying drawings, FIGS. 1 to 3 are graphs showing the results of the measurements conducted by the method and the sensor of the present invention.

FIG. 4 is a graph showing the results of measurement conducted by a conventional method.

FIG. 5 is a block, diagram of a conventional sensor.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

Now, the present invention will be described in detail with reference to the preferred embodiments.

In the first and third embodiments of the present invention, the exponential function of the correction factor is determined by the following procedure.

For instance, 1 km of an optical fiber to be measured is divided into 100 measuring points with an interval of every 10 m. At the first measuring point x=1, the temperature T(1) is calculated by the equation containing no correction factor. The temperature T(2), at the second measuring point x=2, can be obtained approximately by introducing the attenuation difference α{T(1)} obtained from the temperature measured at the first measuring point into the correction factor, since the error is negligible. Further, the temperature T(3), at the third measuring point x=3, can likewise be obtained by introducing the attenuation difference α{T(1)} at the first measuring point and the attenuation difference α{T(2)} at the second measuring point into the correction factor. Accordingly, the correction factor C_(n) at the nth measuring point (where n is an integer of at least 1) ##EQU11## In this case, data of attenuation difference α in a temperature range of from about 0° C. to about 100° C. to which the sensor of the present invention will be applied, are memorized or stored by the signal processing unit prior to the measurement of the temperature distribution and will be used for the measurement. Here, the attenuation difference α is determined by the following procedure. When the temperature of the entire optical fiber is constant, the amplitude signals I_(AS) (x) of the anti-Stokes light is represented by the following equation:

    I.sub.As (x)=A(T)·e.sup.a(T)·x ·(5)

where x is the distance, a(T) is the attenuation of the anti-Stokes light at a temperature T, and A(T) is a function which is dependent on the temperature. When natural logarithm is applied to both sides of the above equation, the following equation (6) will be obtained:

    lnI.sub.AS (x)=lnA(T)+a(T)·x                      (6)

Here, if the temperature T is constant against x, i.e., if the temperature is constant over the entire optical fiber, equation (6) will be a linear line with a slope of a(T) with lnA(T) being intercept. Accordingly, by changing the temperature T to various values and conducting the above treatment, the slopes a(T) at various temperatures can be obtained. The attenuation b(T) of the Stokes light can be obtained in a similar manner. Accordingly, the attenuation difference α(T) at a temperature T will be obtained by the following equation:

    α(T)=b(T)-a(T)(m.sup.-1).                            (7)

Further, as a result of extensive research, the present inventors have found that α(T) can be approximated by a quadric curve relating to the temperature T, and α(T) can be memorized as a quadric function of T, which can in turn be used for the calculation of the correction factor C_(n). Further, α(T) may vary slightly depending upon the type of the optical fiber used, and in such a case, α(T) may be determined for each type of such optical fiber.

Now, the second and forth embodiments of the present invention will be described.

Firstly, the following treatment is conducted by setting u=0 or 1 in the first correction equation. However, the case where u=1 is preferred, since the measurement can thereby be conducted with high precision. In the case of u=1, the attenuation difference α is substituted as a constant value corresponding to a certain temperature within a temperature range of from about 0° C .to about 100° C. in which the sensor of the present invention will be applied. In this case if an optical fiber of 1,000 m is divided into 1,000 measuring points, spaced every 1 m, the temperature distribution [T_(n) ^(o) ] (wherein n is an integer of 1≦n≦1,000) is obtained with respect to the 1,000 measuring points. Also, in the case of u=0, [T_(n) ^(o) ] is obtained in a similar manner. Then, the above temperature distribution is substituted into the exponential part of the second correction equation, i.e., by using the above equation (3), calculation, by the following equation (8): ##EQU12## is conducted with respect to n=1 to 1,000. The result is represented by [T_(n) ¹ ] (1≦n≦1,000). In equation (8), the sum of up to n terms is taken to increase the precision. Next, the error ε is obtained by the following equation: ##EQU13## If the error ε is larger than the prescribed value e.g. 3, [T_(n) ¹ ] is further substituted into the second correction equation to obtain [T_(n) ^(2]) (1≦n≦1,000). This operation will then be repeated until the error ε becomes not higher than the prescribed value, and the temperature distribution [T_(n) ^(r) ], when the error ε becomes not higher than the prescribed value, is taken as the measured value. This operation may be represented by the following equation (10):

    [T.sub.n.sup.n ]=f[T.sub.bn.sup.r-1 ] ;                    (10)

where f[T_(n) ^(r-1) ] is the operation of substituting the previous temperature distribution [T_(n) ^(r-1) ] into the expenential part of the right side of the second correction equation. Namely, the results of the left side is are sequentially substituted into the right side, and the error ε is determined until the desired value is obtained. In the foregoing instance, the measuring points are taken every 1 m, but the measurement may be conducted for every 10 m or other length. Further, in this case, the attenuation difference α of the exponential part of the second correction equation is taken in consideration of the temperature dependency, and its data are memorized or stored by the signal processing unit prior to the measurement of the temperature distribution or they are memorized in the form of a quadric function of the temperature T, and will be used for the measurement of the temperature distribution.

In the present invention, a correction factor, based on the attenuation reference α of the Stokes light and the anti-Stokes light, is introduced into the Dakin equation used for the measurement of the temperature distribution in the optical fiber, and the correction equation taking the temperature dependency of the attenuation difference α into consideration is used, whereby in the temperature distribution measurement of an optical fiber of a long distance, the problem that the temperature distribution curve inclines substantially due to the distance, can be solved, and the measurement can be conducted with high precision.

Now, the present invention will be described in further detail with reference to the following Examples. However, it should be understood that the present invention is by no means restricted to such specific Examples.

EXAMPLES

FIGS. 1 to 3 illustrate the Examples of the present invention.

FIG. 1(a) is a graph showing the temperature distribution obtained by using the correction equation: ##EQU14## FIG. 1(b) is a graph of a Reference Example wherein α is a constant value, and the exponential part is simply represented by αx. In FIG. 1(a), the temperature distribution is not substantially inclined as compared with the conventional instance of FIG. 4, and the precision is substantially improved over the case of FIG. 1(b). In this Example, an optical fiber of 1,000 m was measured at measuring points set every 1 m.

FIG. 2(a) is a graph showing the temperature distribution obtained by the convergent method employing the first correction equation (u=1) having αx introduced into the exponential part and the second correction equation, and it shows the results up to the second treatment, i.e., [T_(n) ² ]. FIG. 2(b) is a graph of a Reference Example wherein α is a constant value, and the exponential part is simply represented by αx. In FIG. 2(a), the temperature distribution is not substantially inclined as compared with the conventional instance of FIG. 4, and the precision thereof is substantially improved over FIG. 2(b). In this Example, an optical fiber of 1,000 m was measured at measuring points of every 1 m. In this Example, errors ε₁ to ε₆ obtained by the operations of 6 times are as follows: ##EQU15## FIG. 2(a) shows a case wherein the operation was terminated when the error became ε₂ ≦3. In the range of ε₃ or below, no difference was visually observed when compared with the graph in the case of ε₂. The prescribed value to be used for the judgment of the error ε is not limited to the above identified value 3 and may suitably be selected.

FIG. 3(a) is a graph showing the temperature distribution obtained by the convergent method employing the first correction equation (u=0) having no exponential part and the second correction equation, and represents the results up to the third treatment i.e., [T_(n) ² ]. FIG. 3(b) is a graph of a Reference Example in which α is a constant value, and the exponential part is represented simply by αx. In FIG. 3(a), the temperature distribution is not substantially inclined as compared with the conventional instance of FIG. 4, and the precision is substantially improved over FIG. 3(b). In this Example, an optical fiber of 1,000 was measured at measuring points of every 1 m. In this Example, the operation was terminated when the error become ε₃ ≦3,whereby the precision was slightly lower than the case of FIG. 2(a). The slightly lower precision is believed to be attributable to the fact that no exponential part was introduced into the first correction equation.

Here, the sensor of the present invention is such that in the conventional apparatus as shown in FIG. 5, the treatment having the above correction factor introduced is conducted by a signal processing unit i.e. a unit comprising a filter 13, photo-electric converters 14 and 14', preamplifiers 15 and 15', an averager 16 and a signal processing section 17. More specifically, the treatment having the above correction factor introduced, is conducted by the signal processing section 17 or by a microcomputer connected to the signal processing section 17.

The results of FIGS. 1 to 3 are those obtained by the measurements wherein α was 0.6×10³¹ 5 (m³¹ 1, Θ was 25(° C.), R'(T) and R'(Θ) were empirically obtained, k=1.38065800×10⁻²³ J.K⁻¹, h=6.62607550×10⁻³⁴ J·S, c=2.99792458×10⁸ m·s⁻¹, and ν=440 (cm¹).

Further, ν may have a varied value or may be an empirically obtained real value. The adverse effect to the relative ratio of the Stokes light and the anti-Stokes light caused by the optical connection between the respective units in the sensor system of the present invention, may be corrected by subjecting the signal processing unit to off setting adjustment by means of a new correction factor δ as shown below, prior to the measurement: ##EQU16##

As described in the foregoing, according to the present invention, treatment is conducted based on an equation having a correction factor introduced to the conventional Dakin equation and further a correction equation taking the temperature dependency of the attenuation difference α in the correction factor into consideration, is employed, whereby the conventional problem such that the entire curve of the temperature distribution substantially inclines with distance, can be solved, and the precision in the temperature measurement for a long distance can be improved. 

We claim:
 1. A temperature-measuring method which comprises the steps of:maintaining a reference temperature point of an optical fiber at a known reference temperature; inputting laser pulses into an incident end of the optical fiber to be measured; detecting the amplitude and delay time of the Stokes and anti-Strokes light contained in a return beam from the optical fiber; and measuring the temperature distribution in the optical fiber from a ratio of respective amplitudes of said Stokes and anti-Stokes light and said delay time of the Stokes light and the anti-Stokes light contained in said return beam from the optical fiber, wherein the temperature distribution is measured by using the following first and second correction equations: first correction equation: ##EQU17## where T(x) is the temperature to be measured at a measuring point, Θ is the reference temperature at the reference temperature point, R'(T) is the relative ratio of the amplitudes of the Stoke and anti-Stoke light of the return beam at the reference temperature point, k is the Boltzmann's constant, h is the Planck's constant, c is the velocity of light, ν is the Raman shift, α is an attenuation difference in the optical fiber between the Stokes light and the anti-Stokes light, x is a distance of the measuring point from the incident end of the optical fiber, and u is a real number of at least 0, wherein the attenuation difference α is represented by a constant value, and an exponential function of the first correction equation in which the exponential part is the product of the attenuation difference α and the distance x, is employed as a correction factor; and second correction equation: ##EQU18## where T(x), Θ, R'(T), R'(Θ), k, h, c, ν, α and x are as defined above, wherein the attenuation difference α is represented by a function α{(T(τ)} which is dependent on the temperature (t(τ) at the measuring point τ, and an exponential function of the second correction equation in which the exponential part is represented by a valued obtained by integrating the function α{T(τ)} with a variation dτ in distance, is employed as a correction factor; in such a manner that at first, temperature distribution T_(n) ^(o) (wherein n is an integer of at least 1) obtained by the first correction equation, is substituted into the exponential part of the second correction equation, to obtain temperature distribution T_(n) ¹, and subsequently, temperature distribution T_(n) ^(r) (wherein r is an integer of at least 1) obtained immediately before is successively substituted into the second correction equation to obtain temperature distribution T_(n) ^(r+1) at r+1 times, and temperature distribution T_(n) ^(r) obtained when an error ##EQU19## (wherein m is the number of measuring points), becomes not higher than a predetermined value at r times, is taken as a measured value.
 2. The temperature-measuring method according to claim 1, wherein the error ε is more than
 3. 3. The temperature-measuring method according to claim 1, wherein prior to the measurement of the temperature distribution, α at various temperatures is memorized by the signal processing unit.
 4. A distributed optical fiber temperature sensor which comprises: means for maintaining a temperature reference point of an optical fiber at a known reference temperature, a light source which inputs laser pulses into an incident end of the optical fiber to be measured, and an optical directional coupler which introduces a return beam from the optical fiber to a signal processing unit, the signal processing unit including means for detecting the Stokes light and the anti-Stokes light contained in the return beam and means for measuring the temperature distribution in the optical fiber from the ratio of the amplitudes and the delay time of the Stokes light and the anti-Stokes light, wherein said temperature distribution is obtained by suing the following first and second correction equations:first correction equation: ##EQU20## where T(x) is the temperature to be measured at a measuring point, θ is the reference temperature at the reference temperature point, R'(T) is the relative ratio of amplitudes of the Stokes and anti-Stokes light at the measuring point, R' (Θ) is the relative ratio of amplitudes of the Stokes and anti-Stokes light at the reference temperature point, k is the Boltzmann's constant, h is the Planck's constant, c is the velocity of light, ν is the Raman shift, α is an attenuation difference in the optical fiber between the Stokes light and the anti-Stokes light, x is a distance to the measuring point from the incident end of the optical fiber, and u is a real number of at least 0, wherein the attenuation difference α is represented by a constant value, and an exponential function of the first equation in which the exponential part is the product of the attenuation difference α and the distance x, is employed as a correction factor; and second correction equation: ##EQU21## where T(x), Θ, R'(T), R'(Θ), k, h, c, ν, α and x are as defined above, wherein the attenuation difference α is represented by a function α{T(τ)} which is dependent on the temperature T(τ) at the measuring point τ, and an exponential function of the second correction equation in which the exponential part is represented by a valued obtained by integrating the function α{T(τ)} with a variation dτ in distance, is employed as a correction factor; in such a manner that at first, temperature distribution T_(n) ^(o) (wherein n is an integer of at least 1) obtained by the first correction equation, is substituted into the exponential part of the second correction equation, to obtain temperature distribution T_(n) ¹, and subsequently, temperature distribution T_(n) ^(r) (wherein r is an integer of at least 1) obtained immediately before is successively substituted into the second correction equation to obtain temperature distribution T_(n) ^(r-1) at r+1 times, and temperature distribution T_(n) ^(r) obtained when an error ##EQU22## (wherein m is the number of measuring points), becomes not higher than a predetermined value at r times, is taken as a measured value.
 5. The distributed optical fiber temperature sensor according to claim 4, wherein the error ε is not more than
 3. 6. The distributed optical fiber temperature sensor according to claim 4, wherein prior to the measurement of the temperature distribution, α at various temperatures is memorized by the signal processing unit. 